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SWITCHING  SYSTEMS,  I1 


Thomas  I.  Seidman2,  3 


ABSTRACT: 

A  model  is  introduced  for  a  switching  system  consisting  of  a  number  of  modes 
(e.g.,  differential  equations)  together  with  a  set  of  switching  rules.  The  rules  per¬ 
mit  occasional  nonuniqueness  but  are  shown  to  preserve  a  continuous  dependence 
property:  the  limit  of  solutions  is  a  solution. 
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1.  Introduction 

By  a  switching  system  we  mean  a  (finite)  set  of  modes  together  with  a  set  of 
switching  rules  of  a  special  form.  A  more  detailed  definition  will  be  provided  later 
but,  for  orientation,  let  us  consider  a  prototypical  example. 

Example  1:  Let  X  be  the  plane  SR2  and  suppose  one  has  available  two  modes  given 
by  the  differential  equations 

x  =  fj{x)  (j  =  1,2;  xEX) 

(Actually,  we  suppose  that  there  are  (closed)  sets  Uj  C  X  for  which  each  mode  is 
available  so  it  is  only  in  Ui  D  U2  that  both  j  =  1  and  j  =  2  would  be  available.)  The 
switching  rule  is  a  modified  form  of  the  heuristic:  Don’t  switch  unless  you  must!  In¬ 
troducing  the  complementary  forbidden  sets  Rj  :=  X\Uj,  we  consider  the  interesting 
case  in  which  Ri  contains  a  global  attractor  for  the  mode  x  —  fi{x)  and  similarly  for 
R 2',  we  assume,  here,  that  R i,  R2  are  disjoint  so  at  least  one  mode  is  always  available. 
Suppose  we  start  the  system  with  state  x(0)  in  (the  interior  of)  U\  n  U2  and  in  mode 
1.  Our  solution  coincides  with  the  solution  of  x  =  fi(x)  until  this  hits  the  boundary 
dlli  =  dRi]  this  is  inevitable  since  we  have  assumed  £1  contains  a  global  attractor. 
Since  it  is  forbidden  to  enter  Ri  in  mode  1  (i.e.,  for  a  state  x{t)  G  Ri  the  mode  j  =  1 
is  unavailable),  we  switch  to  mode  2.  With  this  change  in  mode  the  trajectory  may  or 
may  not  enter  £1  but  coincides  with  a  solution  of  x  ~  f2(x )  until  hitting  dU2  —  dR2. 
We  expect,  then,  to  follow  the  individual  differential  equations  alternately,  switching 
from  j  =  1  to  j  =  2  at  dRi  and  back  from  j  =  2  to  j  =  1  at  dR2 ■  □ 

An  apparently  minor  quibble  arises:  Suppose  there  would  be  a  point  £  G  dRi  for 
which  the  equation  x  =  fi{x)  has  a  trajectory  t  passing  through  £  but  staying  locally 
in  Ui  (see  Figure  1),  i.e., 
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tangential4  to  dR  i  at  £.  Do  we  or  do  we  not  switch?  This  figure  shows  a  neighboring 
trajectory  t  entering  R\  at  £1  so,  for  the  switching  system,  the  policy  is  clear:  switch 
modes  at  £j.  On  the  other  hand,  the  figure  shows  another  neighboring  trajectory  r2 
which  never  hits  dR\  until  £2  so  the  policy  for  the  switching  system  is  again  clear: 
switch  modes  at  £2.  Taking  limits  through  trajectories  like  T\  we  would  expect  switch¬ 
ing  at  £  for  a  switching  system  trajectory  coinciding  (in  part)  with  the  earlier  part 
of  t  while  taking  limits  through  trajectories  like  r2  would  suggest,  for  the  switching 
system,  that  one  continue  to  coincide  with  r  and  defer  switching  until  £'.  In  order 
to  have  any  chance  of  preserving  the  principle  that  a  limit  of  solutions  should  be  a 
solution,  we  accept  both  possibilities:  switching  at  £  and  deferring  switching  until 
£'.  This  means,  of  course,  that  we  must  accept  the  consequence  that  a  solution  for 
the  switching  system  which  initially  coincides  with  the  early  part  of  r  must  then 
have  a  nonunique  continuation.  This  possibility  of  nonuniqueness  is  a  significant 
characteristic5  of  the  theory  of  switching  systems  —  although  for  the  considerations 
of  this  paper  it  does  not  materially  affect  the  results.  Even  apart  from  this  possibility, 
we  will  see  from  Examples  3,  and  4  that  care  is  needed  in  deducing  properties  of  a 
switching  system  from  corresponding  properties  of  the  modes  used  to  define  it. 

In  the  next  section  we  will  (very  briefly)  indicate  some  possible  applications/ 
examples  of  the  formulation  although  a  more  detailed  exploration  of  these  will  be 
deferred  to  a  later  presentation.  Here,  these  are  only  intended  as  suggestive  and  to 
motivate  the  precise  formulation  of  the  switching  rules  in  Section  3. 

Section  4  contains  our  principal  result,  on  the  continuous  dependence  of  regular 
solutions,  together  with  a  global  existence  result.  Section  5  discusses  an  important 
special  case:  linear  switching  systems.  Finally,  the  last  section  will  note  some  open 
problems  and  directions  for  generalization.  A  principal  topic  for  further  discussion  is 
the  existence  of  periodic  solutions  but  this  will  be  deferred  to  a  separate  paper  [6]; 
it  is  in  this  context  that  the  possibility  of  nonuniqueness  described  above  becomes 
overwhelmingly  significant. 

2.  Motivation 

The  original  motivation  for  formulating  a  notion  of  ‘switching  systems’  came  from 
an  attempt  to  model  thermostats.  The  two  modes,  in  this  case,  correspond  to  FUR¬ 
NACE  OFF  and  FURNACE  ON  and  the  thermostat  is  a  device  to  switch  between  these. 

4  Although  shown  as  such  in  Figure  1,  this  need  not  mean  tangency  in  the  usual  geometric  sense  since 
we  do  not  impose  enough  regularity  on  dZj  for  this  to  be  necessarily  meaningful;  it  might  be  better 
to  say  that  r  fails  to  be  transversal  at  £.  All  we  really  mean,  here,  in  using  the  term  ‘tangential’,  is 
that  r  contains  a  point  f  of  dZ i,  but  does  not  enter  Z\  at  £. 

6Note  that  the  geometry  of  Figure  1  is  generic  in  that  some  (nearby)  point  of  tangency  must  occur 
for  any  small  perturbation  of  the  direction  field  and/or  the  boundary  8Z\. 
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It  consists  of  a  sensor  (measuring  the  temperature  9  at  a  particular  position)  and  a 
pair  of  set  points  9i,  82  (Typically,  the  gap  02  —  #i  >  0  is  fixed  and  the  mean  ($i  +  $2) / 2 
is  adjustable.)  with  two  internal  states  corresponding  to  the  two  modes.  If  the  fur¬ 
nace  is  OFF,  then  it  will  be  switched  ON  when  9  crosses  9\  from  above;  in  particular, 
the  furnace  will  always  be  on  when  9  <  9 1.  When,  eventually,  the  temperature  rises 
to  have  6  cross  the  upper  set  point  62  (from  below),  then  the  furnace  will  be  switched 
OFF;  no  recrossing  of  the  lower  set  point  affects  the  state  until  this  occurs;  the  furnace 
is  always  OFF  for  0  >  02. 

This,  together  with  the  partial  differential  equation  governing  the  evolving  state 
(spatial  temperature  distribution),  seems  to  describe  the  physics  quite  well  except  for 
the  same  minor  quibble  noted  in  the  Introduction:  What  happens  if,  e.g.,  with  the 
furnace  ON,  the  sensed  temperature  9{t)  rises  to  62  without  (immediately)  crossing? 
(Since  the  evolution  is  given  by  a  pde,  one  can  find  initial  conditions  for  which  this  is 
would  actually  occur.)  It  was  the  analysis  of  this  situation  which  led  to  the  present 
model.  A  completely  different  model  is  discussed,  for  example,  in  [4];  the  discussion 
in  [1]  is  closer  in  its  concerns  with  the  present  analysis.  More  detailed  analysis  of 
switching  system  thermostat  models  will  be  deferred  to  [7]. 

An  entirely  different  setting  leading  to  switching  system  models  is  the  Hamilton 
-  Jacobi  -  Bellman  formulation  of  optimal  feedback  control  of  multimodal  (variable 
structure)  systems  with  switching  costs.  Suppose  one  has  a  system  which  can  operate 
in  any  of  J  modes,  e.g.,  corresponding  to  differential  equations 

(2.1)  x  =  fj(x)  = 

with  attendant  (x,  independent  running  costs  but  with  the  possibility  of  switching 
at  any  time  from  the  i~th  to  the  k-th  mode  with  cost  Cjik(x)\  there  may  also  be  other 
control  possibilities  implicit  in  (  2.1).  We  suppose  Vj(x)  is  the  optimal  infinite  horizon 
(discounted)  cost  if  one  is  at  the  state  x  in  the  mode  j.  Clearly,  we  would  switch  to 
the  A>th  mode  if  V)(x)  >  Vk(x)  +  Cjtk(x)  but  would  not  switch  if  the  reverse  inequality 
would  hold,  i.e.,  we  would  always  choose 

(2.2)  k  =  argmin{Vk(x)  +  cJ|fc(x)  :  k  — 

where  ( x,j )  is  the  current  state  and  k  is  to  be  the  ‘new’  mode  (set  Cjj(x)  :=  0  for 
completeness).  For  a  discussion  of  this  approach,  see  [3].  The  effort,  then  has  gone 
into  the  construction  of  the  value  functions  {V)(-)}  with  (  2.1,  2.2)  taken  as  defining 
the  controlled  dynamics.  We  observe  that  the  possibility  of  nonuniqueness  in  (  2.2) 
means  that  a  further,  more  detailed,  analysis  of  the  dynamics  is  needed  precisely  at 
the  switching  surfaces 


Sj,k  :=  {x  :  Vj(x)  =  Vk(x)  +  cjik(x)}. 
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If  we  let  Rj  be  the  open  set 

Zj  :=  {x  :  Vj(x)  >  Vk{x)  +  cjtk(x)  for  some  k  ^  j}, 

then  we  will  obtain  a  switching  system  model  for  the  optimally  controlled  dynamics. 
Note  that  the  possibility  of  nonunique  continuation  if  a  trajectory  r  of  x  —  /y(x)  is 
tangential  to  the  switching  surface  simply  means  that  the  optimal  cost  is  attained 
nonuniquely:  either  continuing  without  switching  or  paying  the  switching  cost  and 
continuing  in  a  new  mode  give  the  same  (optimal)  cost. 

There  appear  to  be  connections  between  the  theory  of  switching  systems  presented 
here  and  ‘viability  theory’,  c.f.,  [2].  Of  particular  interest,  in  this  connection,  is  the 
notion  of  a  ‘heavy  trajectory’,  cf.  [  ],  corresponding  to  the  reluctance  to  change  modes 
implicit  in  our  switching  rules.  Our  considerations  are,  however,  almost  disjoint  from 
those  of  [2]  since  we  can  also  write,  e.g.,  (  2.1)  as  a  differential  inclusion  x  6  F(x) 
by  setting  F(x)  :=  {fj(x)  :  j  —  1, . . . ,  J}  but  we  are  obviously  emphasizing  the  case: 
F(x)  finite  whereas  [2]  emphasizes  the  quite  distinct  case:  F(x)  convex. 

A  final  motivating  setting  comes  from  singular  perturbation  theory.  Consider,  for 
example,  a  system 

(2.3)  x  =  /(x,y),  ey  =  g{x,y) 

for  very  small  e:  >  0.  A  principal  concern  of  singular  perturbation  theory  is  to  initiate 
an  analysis  of  (  2.3)  by  comparison  with  the  reduced  order  (implicit)  model 

(2.4)  x  —  f(x,y)  ,  y(x,y)  =  0. 

If  we  can  solve  g(x,y)  =  0  to  obtain  y  =  Y(x),  then  (  2.5)  becomes  simply 

(2.5)  x  =  f*{x)  :=  f{x,Y{x)). 

Suppose,  however,  the  graph  {(x,y)  :  g{x ,  y)  =  0}  would  look  like  that  of  Figure  2. 
As  shown,  we  observe 


Figure  2 
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that  A  =  dg/dy  is  negative,  corresponding  to  stability  of  the  perturbation  equation: 
ey  =  Ay,  along  the  branches  y  =  Yj(x )  for  j  =  1,2  but  we  have  A  >  0  (instability) 
on  y  —  Y*(x).  Thus,  local  analysis  shows  that  once  we  have  y  ps  Yi(x )  with  a  <  x  we 
would  expect  to  stay  close  to  the  solution  of 

x  =  fi(x)  :=  f{x,Y i(x)) 

(for  very  small  e  >  0)  unless/until  this  solution  would  reach  x  =  a.  Assuming 
/i(«)  <  0,  this  trajectory  would  (try  to)  enter  the  forbidden  region  R i  :=  (— oo,a)  for 
which  Yi  does  not  exist.  One  would  have  to  make  a  more  detailed  analysis  of  (  2.3) 
to  see  just  what  would  happen  then  but,  assuming  Y2(x)  would  be  globally  attractive 
for  <f  =  g(x,  f)  for  x  <  a,  we  would  expect  a  rapid  transient  behavior  ‘switching’  to 
the  branch  y  =  y2(x).  The  solution  would  then  stay  close  to  this  branch  unless/until 
it  reached  x  =  b,  with  a  similar  ‘switching’  there  if  /2 (b)  >  0.  This  behavior  is,  of 
course,  the  standard  hysteresis  loop  and  is  essentially  similar  to  the  description  of 
the  thermostat  above.  Other  interesting  varieties  of  behavior  become  possible  when 
x,y  may  ‘live’  in  higher-dimensional  spaces.  We  defer  to  [6]  any  further  discussion  of 
this  and  of  the  extent  to  which  the  switching  system  model  may  provide  a  suitable 
reduced  order  model  with  which  to  initiate  the  analysis  of  (  2.3).  The  considerations 
seem  complementary  to  those  of  [5] . 

3.  Definitions 

Motivated  by  the  examples  of  the  previous  section,  we  introduce  a  formal  defini¬ 
tion  of  ‘switching  system’.  This  may  not  be  the  most  general  possible  notion  (c.f., 
Section  6  for  further  comment)  but  more  than  covers  the  motivating  settings. 

First,  by  a  mode  we  mean  a  semidynamical  system  with  state  space  X: 

(3.1)  7r  :  5R+  x  X  — >  X  continuous,  with  7r(0,  £)  =  £  and 

7r (t  +  s,  £)  =  7r(f ,  tt(.s,  £))  for  t,  s  >  0. 

Note  that  tt  may  be  interpreted  as  the  solution  map:  £  =:  x(0)  i — *  x(t)  for  an 
autonomous  differential  equation6  x  =  f(x)  but  the  exposition  is  somewhat  simpler 
this  way  and  more  general  in  some  respects.  We  will  consider  a  variable  structure 
system  corresponding  to  a  set  {7 r;-  :  j  G  J}  of  such  modes  where  J  is  suitable  index 
set  (usually  finite:  j  =  1, . . . ,  J  ). 

For  each  j  £  J  we  assume  we  are  given  a  forbidden  set  Rj  C  X  while  for  each  pair 
(j,  k)  e  J  x  J  there  is  an  admissible  switching  set  Sj,k  C  X.  (Any  of  these  sets  may 

6This  restriction  to  autonomous  equations  is  not  very  significant  in  the  sense  that  one  can  formally 
include  t  as  a  component  of  the  state.  Somewhat  more  restrictive  is  the  implicit  assumption  that  a 
unique  global  solution  (i.e.,  for  all  t  6  5R+)  exists  for  each  initial  £  €  X,  see  Section  6. 
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possibly  be  empty.)  Our  concern  with  Rj  is  that  the  mode  7 ry  may  possibly  enter  Rj 
at  some  £  £  dRj,  by  which  we  mean 

(i)  for  some  £0  we  have  7 ry(-,  £0)  6  X  \  Rj  on  (0,  er)  with  e  >  0 

(3.2)  and  7T,-(e,  £0)  =  £\ 

(ii)  7 Tj(e,  £)  G  for  arbitrarily  small  e  >  0. 

We  impose  the  set  of  geometric  hypotheses:7 

0)  $j,k  0  Rk  =  0  ,  Rj  =  X  \  Sj.j ; 

(ii)  each  Sjtk  is  closed  in  X  (so  each  Rj  is  open); 

(3.3)  (iii)  for  each  £  G  dRj  there  is  some  k  G  J  such  that  £  G  Sj,k] 

(iv)  for  each  j  G  J  :  for  each  £  G  X  there  is  a  neighborhood  .A/y ( £) 
which  nontrivially  intersects  only  finitely  many  of  { Sj ,&}. 

By  a  switching  system  we  mean  a  specification: 

S  :=  [J,  {7 Tj  :  j  G  J},{Sjj k  :  /,  A;  G  J}] 

subject  to  (  3.1)  and  (  3.3). 

For  a  function:  t  t— >•  [x(t),j(t)\  :  9i+  — >  X  x  J  to  be  considered  as  a  possible 
solution  of  such  a  switching  system  E  we  first  require  that 

(i)  j(-)  is  piecewise  constant  with  isolated8  ‘jumps’; 

(3.4)  (ii)  whenever  j(-)  is  constant  (=  k)  on  an  interval  (s,t)  we  have 

z(£)  =  Kk{t  -  s,  x(s)). 

Note  that  this  automatically  forces  x(-)  to  be  continuous:  9?+  — ►  X  since  each  mode 
7 Tk  is  continuous  in  t.  The  requirement  that  the  ‘jumps’  in  j (■)  be  isolated  means  that 
jf'(-)  will  have  both  left-  and  right-handed  ‘limits’  at  a  jump: 

(3.5)  j{t-)  =  j  ,j{t+)  =  k  (with  j  ^  k). 

We  refer  to  such  a  time  t  as  a  switching  time  and  to  the  situation  as  a  transition 
(from  mode  iXj  to  mode  7 r*;)  or,  more  succinctly,  as  a  switch:  j  k.  If  t'  is  a  switching 
time  (and  there  are  any  subsequent  switches)  then  there  is  a  unique  next  switching 

7  This  set  might  be  weakened  slightly  if  we  did  not  prefer  to  keep  it  independent  of  any  knowledge 
of  { 7Ty } .  For  example,  (iii)  is  really  needed  only  if  7 ry  would  enter  Rj  at  £  and  even  then  we  might 
permit  £  E  dRk  if  7 a  does  not  enter  Rk  at  £;  compose  the  reformulation  in  Section  6. 

8This  requirement  could  be  relaxed  somewhat  and  we  could  refer  to  the  notion  here  as  defining  a 
regular  solution  of  E. 
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time  t".  We  refer  to  as  an  inter  switching  interval  9  We  impose  the  switching 

rules : 


(i)  a  switch:  j  ^  k  is  permitted  at  time  t  only  if  x(t )  6  Sy,* 

(ii)  if  x(t)  G  £y,  then  j(t)  =  j  is  forbidden  —  hence  t  must  be  a 
switching  time  if  7 ry  enters  £y  at  x(t)  when  one  ‘arrives’ 

at  x(t)  =  £  G  3£y  with  j(i— )  =  j. 


Thus,  by  a  solution  of  the  switching  system  E  on  the  time  interval  [0,T)  we  mean 
a  function  pair 

MO,  j(-)| :  |0,T)  — >Xx  J 


subject  to  (  3.4)  and  (  3.6).  For  definiteness  we  take  j(-)  to  be  left-continuous.  The 
values  [:r(0),  j(0)]  are  called  the  initial  data  for  the  solution;  we  have  x  continuous  at 
0  so  x(0+)  =  a:(0)  but  we  do  permit  (subject  to  (  3.6)  (i))  an  immediate  switch  for 

/(•). 


Remark  1:  Usually  we  only  consider  solutions  with  T  =  oo,  i.e.,  [0,  T)  :=  3J+.  It 
should  be  clear  that  there  is  nothing  ‘magic’  about  0  and  any  interval  I  :=  [a,  b )  would 
be  equally  appropriate  but,  as  we  consider  only  autonomous  systems  here,  there  is  no 
loss  of  generality  in  translating  in  time  by  a  to  consider  [0,  T)  with  T  :=  b  —  o;  initial 
data  at  t  =  a  then  becomes  initial  data  at  t  =  0.  The  other  interesting  possibility 
would  be  to  consider  solutions  ‘for  all  time’  :  1  =  3?;  in  this  case  it  is  somewhat 
irrelevant  to  consider  ‘initial  data’.  This  possibility  will  be  of  interest  to  us  in  the 
context  of  stable  linear  switching  systems  in  Section  5:  note  Theorem  6. 


Example  2:  We  present  here,  partly  for  its  own  sake10  and  partly  in  contrast  with 
the  flavor  of  Example  1,  an  example  with  an  infinite  dimensional  state  space  X. 

We  take  this  space  to  be11 


(3.7)  :={£(•)  :  measurable  from  3R+  to  [0,1]}. 

The  metric  topology  we  consider  on  X°  will  be  that  induced  by  a  weighted  jLx-norm: 

r  oo 

(3.8)  <*(£,£')  :=  /  £>(s)|f(s)  -  £'(«)[ ds, 

Jo 

9We  have  j(  )  constant  (say,  =  k)  on  such  an  interswitching  interval  and  (  3.4)  (ii)  requires  the 
solution  simply  to  follow  the  mode  Wk  during  it. 

10The  specific  significance  of  this  particular  example  lies  in  its  canonical  relation  to  general  linear 
switching  systems  as  discussed  later  in  Theorem  6. 

11  The  particularly  interesting  subset  of  X°  is 

X0  ■—  {£(■)  :  measurable  from  5R+  to  {0,  l}} 
and  we  could  really  have  used  any  X  suitably  embedding  XQ ■ 
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where  we  assume  that: 


(3.9) 


the  weight  (p  is  positive  and,  nonincreasing  with 


e.g.,  <p{s)  =  Ce  as. 
(3.10) 


r  oo 

/  (p(s)ds  =:  M  <  oo, 
Jo 


Next,  set  J  :=  {1,2}  and  define  the  modes  TTj (j 

for  s  >  t, 

j  -  1  on  [0,  t] 


b(«.e)]W=='f  f(s  () 


1,2)  by 


for  t  £  5J+  and  £  £  X°.  One  easily  verifies  that  each  n j  satisfies  (  3.1).  Even  before 
introducing  the  sets  {S3lk}  we  note  a  representation  formula 


(3.11) 


[*(*)](«) 


£(s  —  t)  for  s  >  t, 

j(t  -  s)  -  1  on  [0 ,f] 


for  any  measurable  switching  function  j(-),  not  necessarily  satisfying  (  3.4)  (i).  Note 
that  this  automatically  makes  [x(-),  ./(•)]  satisfy  (  3.4)  (ii)  whenever  (  3.4)  (i)  really 
does  hold.12 

Next,  suppose  we  have  a  specified  function  ip  :  — >  8?  such  that: 


(i)  V>(0)  =  0,  /0°°  i>(s)ds  =  1; 

(3.12)  (ii)  sup{|^(s)|  :  s  >  t}  <  <p{t)-, 

(iii)  \ip[t)  —  0(a)|  <  K<p(s)(t  —  s)  for  t  >  s  >  0. 


We  then  introduce  the  sensor  functional 

r  oo 

(3.13)  £  ►  9{£]  ■=  /  0(b)  Z{s)ds 

Jo 

and,  for  any  X°-valued  state  function  a:(-),  the  associated  sensor  function 

*  ^  r  oo 

(3.14)  t  i — *•  $(t)  =  9{t\x)  0[x(t)]  :=  /  ^(«) [a:(i)] (s)  ds. 

Jo 

Note  that  9  is  continuous  on  X°  and  that  if  x(-)  is  given  by  (  3.11)  we  have 

A  f  OO  /*t 

(3.15)  6(t)  =  /  ip(t  +  s)£(s)ds  +  /  ip(t  —  s)[j(s)  —  l]ds. 

Jo  Jo 

so  that,  by  (  3.9)  and  (  3.12)  (ii),  6{t)  is  almost  independent  of  £  for  large  t :  the 
^-dependent  part  of  9{t)  is  the  first  integral  in  (  3.15)  which  is  bounded  by  /t°°  0 


12Observe  also  that  Xo  is  invariant:  if  £  =  2(0)  €  Xo,  then  x(t)  £  Xo  for  each  t  £  5R+  for  any  j{  ). 
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which  goes  to  0  as  t  —*■  oo.  We  have 
a  a  r°° 

6(t)—6(t')  =  /  {ip(t  +  s)  —  +  s)]£(s)ds 

Jo 

+  [  [ilf{t-s)-il)(t' -s)]\j(a) -l]ds 
Jo 

+ L  ^  -  i^s 

for  t  >  t'  >  0  so,  using  (  3.12)  (ii,  iii),  we  obtain  the  uniform  Lipschitz  condition: 
a  a  r° o 

(3.16)  |0(f)-0(f')|  <  /  K<p{t' +  s)(t  -  t')ds 

Jo 

rt '  ft—t' 

+  /  K<p(t'  —  s)(t  —  t')ds  +  /  ip(s)ds 
Jo  Jo 

<  [KM +  <p{0)}{t-t'). 

In  terms  of  $[•]  we  define,  finally, 

Zi  =  {£  €  X  :  0[t]  <  Cl}t  Z2:={teX:  0[£]  >  C2}, 

(3.17)  Si,i  :=  {£  G  X  :  0[£]  >  Cx},  S2)2  :=  {£  e  X  :  0[£]  <  C2}, 

Sii2  :=  dR1  =  {£eX:  0[£]  =  CJ,  S2,2  :=  dR2  =  {£  e  X  :  0[£]  =  C2} 

where  Ci,  C2  are  specified  constants  with  0  <  Ci  <  C2  <  1.  One  immediately  verifies 
(  3.3)  since  <?[•]  is  continuous.  Note  that  for  a  (possible)  solution  [a:(-) ,  j  (•)]  we  can 
have  a  switch:  j  a*  k  at  time  t  only  if  6(t)  =  Cj.  Thus,  if  (i„,f„+i)  is  an  interswitching 
interval  we  have  |0(i„+1)  —  <?(£„) |  =  C2  —  C\  so,  by  (  3.16), 

(3.18)  tu+ 1  -tv>  [C2  -  CO /[KM  +  £>(0)]  =:  Tmin 
which  certainly  ensures  (  3.4)  (i). 

It  is  interesting  to  note  that  if,  e.g.,  we  were  to  proceed  in  the  mode  7r2  then 

x(t  +  t)  =  7r2(f,  £)  with  £:=x(r), 


9{r  +  t)  =  <9[?r2(i,  ^)] 

=  f  ip(t  +  $)£(s)ds  +  f  ip(t  —  s)[2  —  l]ds 

Jo  Jo 

>  [  i/j(s)ds  — »  1  by  (  3.12)  (i) 

Jo 

so  R2  is  a  global  attractor  for  7t2  and,  indeed,  one  reaches  R2  (from  any  starting  point 
£  G  X°,  proceeding  by  tt2  in  time  not  greater  than  r  —  t2  where 
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Similarly 

poo 

=  /  ip{t  +  s)Z(s)ds 

Jo 

<  f  ,ip(s)ds  =  1  —  f  ip(s)ds  — >  0 

Jt  Jo 

so  is  a  global  attractor  for  7Ti  and  one  has  an  upper  bound  T\  for  the  time  to  reach 
Si. 

For  the  switching  system  E  :=  [J",  {7Ty},  {^y.*}],  this  ensures  that,  for  every  con¬ 
sistent  starting  point  [z(0),y(0)]  (i.e.,  i(0)  %j(o)  )  one  has  at  least  one  solution  of 

E  —  e.g.,  corresponding  to  switching  whenever  this  is  permitted  by  (  3.6)  (i)  —  and 
every  solution  switches  infinitely  often  with  the  uniform  bounds 

Tynin  ^  ^i/+l  ^1/  5:  Tmax  !=  171Q,x{ti,  T2} 

for  interswitching  intervals.  Clearly,  the  state  function  x(-)  is  always  related  to  the 
switching  function  j(-)  through  the  representation  formula  (  3.11)  and  one  relates  j(-) 
to  x(-)  through  0(-)  and  (  3.6);  one  always  has  (  3.4)  in  this  case.  □ 

We  next  consider  the  appropriate  notion  of  convergence  for  switching  systems. 
We  say  En  — *  E°  providing 

(i)  Jn,  J°  =  J  (slightly  more  generally,  J°  =  limsupj n  ); 

(ii)  if  £  £  then  £  £  Z”  for  all  n  >  n.y(£); 

(iii)  if  £°  is  the  limit  of  a  subsequence  — ►  £°  with  each 

r,m)  £  then  f°  £  S,x, 

(3.19)  (iv)  (  3.3)  (iv)  holds  uniformly  in  n,  i.e.,  for  each  j,  £  there  are 

a  neighborhood  N  =  and  a  finite  subset  K  =  Ky(£)  C  J 

such  that  S”A  fl  JJ  =  4>  for  k  #  K  and  each  n\ 

(v)  7r "(i,  £)  — »•  7r° (t,  £)  locally  uniformly,  i.e.,  uniformly  on  some 
[o,r]  x  X-neighborhood  of  each  ( t ,  £),  for  each  j. 

(This  is  actually  a  form  of  ‘upper  convergence’  for  switching  systems.  It  would  be 
plausible  to  adjoin  the  condition 

£°  £  Sj  k  only  if  there  is  a  sequence  £"  — *  £°with  £n  6  Spk 

as  a  complement  to  (iii)  above.  Since  we  will  only  seek  to  prove  an  upper  semi¬ 
continuity  result  for  the  solution  set,  the  definition  (  3.19)  is  adequate  for  our  pur¬ 
poses.) 

4.  Basic  Results 

In  this  section  we  formulate  and  prove  our  basic  general  results.  The  first  two 
results  give  a  kind  of  continuous  dependence  of  the  solution  set  on  the  initial  data 


12 


and  the  system:  we  show  (under  suitable  hypotheses)  that  the  limit  of  solutions  is  a 
solution  when  such  a  limit  exists  and  that  such  limits  always  exist  for  subsequences 
corresponding  to  the  possible  choices  when  nonunique  continuations  are  permitted 
by  the  switching  rules  (  3.6).  We  also  provide  an  existence  result  in  a  somewhat  more 
restricted  setting. 

For  the  first  results  we  must  introduce  a  suitable  notion  of  convergence  for  se¬ 
quences  {[xn,jn]}.  For  the  state  trajectories  {£"(•)}  we  use,  simply,  uniform13  con¬ 
vergence  on  [0,  T]  —  uniform  on  finite  intervals  [0,T']  if  T  —  oo.  For  the  index 
functions  {}"(•)}  the  relevant  notion  is  again  the  natural  one  but  is  a  bit  more  dif¬ 
ficult  to  describe.  It  is  convenient,  here,  to  think  of  an  index  function  ,;(•)  as  being 
specified  by  its  sequence14  of  switching  times 

(4.1)  0  =:  t0  <  ti  <  . . .  <  tp  :=  T 

and  the  assumed  interswitching  indices  {jv u  =  0, 1, ... ,  £>}  where  j (0)  =  j0  and,  for 
i/  =  1,...,P, 

(4.2)  j{t)  =  ju  for  tv- 1  <t  <tu. 

Suppose  we  have  an  index  function  j°(*)  with  isolated  switching  times  {t°v  :  u  = 
0,  ...,P°}  as  in  (  4.1).  If  P°  =  oo  we  say  that  jn{ •)  — >  j°{-)  if  Pn  — ►  oo  =  P°  (e.g., 
Pn  —  oo  for  n  large  )  and,  for  each  p  =  1,2, . . . ,  one  has: 

(i)  tnu^t0v  as  n  —r  oo; 

'  ‘  (ii)  Ju  =  Ju  for  n  >  nv 

If  P°  <  oo  we  ask  instead  that  Pn  >  P°  for  every  (large  enough)  n  and  that  we  have 
(  4.3)  for  each  v  —  1, . . .  ,P°.  (  We  will  also  require  (  4.3)  (ii)  for  v  =  0  although  the 
initial  index  j0  has  no  direct  effect  on  the  state  trajectory  since  the  condition  (  3.4) 
(ii)  is  effectively  applicable  only  for  nontrivial  intervals.) 

Theorem  1:  Suppose  {Sn}  is  a  sequence  of  switching  systems  on  X  converging  to 
a  switching  system  E°  in  the  sense  of  (  3.19)  and  suppose  [ xn,jn ]  is  a  solution  of  En 
on  [0,T]  for  each  n  —  1,2, —  Let  [m°,  j°]  :  3?+  — ►  X  x  J  with  j°(*)  satisfying  (  3.4) 
(i).  Then,  if  [ xn,jn ]  — >  [z°,i°]  as  above,  [x°,j0]  is  a  solution  of  E°  on  [0,T]. 

PROOF  :  This  is  just  the  assertion  that  the  limit  of  solutions  is  a  solution  (modulo 
independent  verification  of  the  regularity  condition  that  switching  times  are  isolated 
for  the  limit).  The  proof  is  quite  straightforward. 

We  have  already  defined  the  switching  times  {t"  :  v  =  0, ...,  vn  \n  —  0,1,...} 
and  now  define  the  corresponding  switching  points  £”  :=  xn(t”).  By  the  uniform 

13  We  will  always  take  X  to  be  a  metric  space  (or  a  topological  vector  space)  so  such  uniformity  is 
meaningful. 

14Note  that  (  3.4)  (i)  prohibits  having  V  =  oo  for  finite  T  but  when  T  =  o owe  can  have  either  u  <  oo 
or  tu  — ♦  oo  as  u  — *  oo  =:  u. 
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convergence,  we  have  £°(-)  continuous  and,  using  (  4.3),  we  see  that  as 

n  — >  oo  for  each  fixed  =  1, We  must  verify  the  switching  rules  (  3.6)  for 
[x°,j°]  and  also  (  3.4)  (ii),  with  (  3.4)  (i)  given  by  assumption. 

Fixing  v,  we  set  j  :=  i°_x(=  for  large  n  by  (  4.3)  )  and  A;  :=  y°(=  i”)  so 
at  t  :=  one  has  a  switch:  j  k  for  j°(-).  Since  £"  G  Sj)*  for  each  n.  =  1,2,... 

and  £"  — >  we  have  G  S?fc  by  (  3.19)  (iii).  Next,  suppose  one  were  to  have 
£  :=  x°(f)  G  R°  for  some  t  which  is  not  a  switching  time  and  had  j°(t)  =  j.  We  would 
then  have  jn(t)  =  j  for  large  n  since  t  must  be  in  the  interior  of  the  interswitching 
interval  (£”_i,i”)  for  each  large  n  with  1/  fixed.  Hence  :=  xn(£)  — »■  f  would  imply 
G  for  large  enough  n  since,  by  (  3.3)  (ii),  R®  is  open.  But  (  3.19)  (ii)  would 
then  give  G  Rf  which  is  a  contradiction.  Thus  (  3.6)  (ii)  holds  at  non-switching 
times.  Since  x°{-)  is  continuous,  G  R°  would  also  give  x°(t)  G  R°  for  nearby 

times  since  R °  is  open;  that  is  impossible.  Thus  (  3.6)  (ii)  always  holds. 

The  verification  of  (  3.4)  (ii)  for  t,s  in  the  interior  of  an  interswitching  interval 
(£°_i,£°)  is  an  immediate  consequence  of  (3.19)  (iv)  since 

x°(t)  =  lim  xn(t)  =  lirmrk(t  —  s,xn(s)) 

=  nk(t  —  s,lim  £n(«s))  ==  7T* (^  —  s,  ar0(s)) 

and  this  extends  to  the  endpoints  by  continuity.  □ 

Theorem  1  shows  that  a  limit,  when  it  exists,  must  be  a  solution.  To  complement 
this,  we  next  show  that,  in  some  sense,  such  a  limit  always  exists.  A  certain  additional 
condition  is  needed  and  we  say  that  a  set  of  index  functions  (or  the  corresponding 
solutions)  is  uniformly  regular  if,  for  each  finite  T'  (we  consider  only  T'  —  T  if  T  <  oo 
),  there  is  a  minimum  length  r  =  t(T')  for  interswitching  intervals  in  [0,  T']  —  more 
formally,  if 

(4.4)  tv  <  T'  =>■  *„_!  +  r(T')  <  ty 

for  each  [tv-\,tv)  associated  with  j(-)  in  the  set. 

Theorem  2:  Suppose  {S'1}  is  a  sequence  of  switching  systems  converging  to  S°  in  the 
sense  of  (  3.19)  and  suppose  [, xn,jn }  is  a  solution  of  En  on  [0,T]  for  each  n—  1,2,  — 
Assume  {}”}  is  uniformly  regular  on  [0,T]  and  x"(0)  =:  ££  — >■  £$  ,jn( 0)  =  j0.  Then 
there  is  a  subsequence  n(m)  such  that  jn(m)]  converges  to  a  solution  [x°,j°]  of 

S°  with  x°(0)  - 

PROOF  :  We  will  proceed  by  repeated  extraction  of  subsequences,  followed  by  a  Can¬ 
torial  diagonal  argument.  Abusing  notation  slightly,  we  continue  to  write  {[xn, J1’"]}, 
etc.,  for  the  subsequences  at  each  stage.  For  exposition  we  treat  only  the  case  T  —  oo 
(so  ‘solution’  means  a  global  solution) ,  leaving  the  trivial  modification  for  T  <  oo  to 
the  reader;  we  also,  similarly,  assume  lim  Pn  =  oo  although  it  may  or  may  not  turn 
out  that  P°  =  oo.  As  above  we  denote  the  switching  points  by  :=  xn(t ”). 
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The  construction  proceeds  by  recursively  (in  u)  obtaining  ft  and  then  ft  ,  This 
will  determine  [x0 (•) ,  j'° (*)] .  The  construction  will  immediately  give  the  convergence 
[xn,jn]  — >•  [x°,i°]  with  (  3.4)  (i)  following  from  the  assumed  uniform  regularity.  Ap¬ 
plication  of  Theorem  (Tl)  then  shows  [x°,j°]  is  a  solution  of  E°  as  asserted. 

We  have  ft  =  ft  and  £”  — »•  by  assumption  and,  of  course,  t"  =  0  =  ft.  It  may  be 
that  ft  =  ft  —  ft  (  infinitely  often)  but  it  is  also  possible  that  ft  is  (infinitely  often) 
a  switching  time  for  jn(-).  In  the  latter  case,  having  a  switch:  j  :=  ft  ^  k{n)  ft 
we  must  have  ££  G  S*k(ny  By  (  3.19)  (iv)  and  the  fact  that  £"  — >■  £q,  this  means 
k(n)  G  K  —  Kj(Co)  f°r  all  but  finitely  many  n.  Since  K  is  finite,  this  means  there 
must  be  at  least  one  index  k  G  K  such  that  k(n)  —  k  for  infinitely  many  n.  Choose 
ft  =  k  and  extract  the  sequence  for  which  k(n)  =  k  =  ft. 

Working  now  with  this  subsequence  (still  denoted  by  [ xn,jn ]  ),  consider  {£”}  which 
necessarily  has  a  convergent  subsequence  (noting  that  we  will  accept  convergence  to 
oo);  we  take  ft  to  be  the  limit  so,  extracting  this  subsequence,  we  now  have  ft  — ►  ft- 
Note  that  the  uniform  regularity  condition  ensures  that 

ft  —  ft  —  ft  >  mm{l,r(l)}  =:  T\  >  0 

so,  in  the  limit,  ft  —  ft  >  t\.  For  ft  <  s'  <  t'  <  ft  we  have  (with  j  —  ft  =  ft  ) 

(4.5)  !"(<)  =  *?(t  -  j;  ,  a)  -  .  eS) 

uniformly  on  [s',t']  by  (  3.4)  (ii),  (  3.1),  and  (  3.19)  (v).  It  follows  that  this  can  be 
taken  to  define  x°(-)  as  a  continuous  function  on  [io^i]  (on  [^o>°°)  if  <?  =  oo  )  with 
:=  x°(ft)  well-defined  if  ft  <  oo.  We  obtain  (  3.6)  (ii)  for  x°(-)  on  exactly 

as  in  the  proof  of  Theorem  1. 

Assuming  ft  <  oo,  it  is  clear  from  (  4.5)  and  (  3.19)  (v)  that  £"  :=  xn(ft)  — > 
x°(t5)  =:  As  before,  since  each  t 1  is  a  switching  time  we  have  Zi  €  %(n)  for  each 

n  with  j  =  ft  =  ft  and  k(n)  —  j2  =  jn(ft+).  As  before  we  can  select  k  G  Kj(ft) 
and  extract  a  subsequence  such  that  each  ft  =  k  =:  ft-  As  before,  we  extract  a 
subsequence  for  which  ft  — *  ft  with  ft  —  ft  >  min{l,T[t\  +  1)}  =:  r2  so  ft  —  ft>  r2. 
As  before  we  have 

(4.6)  xn(t)  =  ft(t  -  ft,  Q)  —  -  tl,  $)  =:  x°(t) 

uniformly  on  any  C  [ft, ft],  defining  x°(-). 

Note  that  the  switching  times  cannot  ‘bunch  up’  since,  for  any  T1  <  oo,  we  must 
have  ft  >  vt[T')  if  ft  <  T1  so  v  <  T'/r(T').  The  same  argument  given  above  for 
construction  of  t]*,  x°(-)  on  [fo,£°],  and  shows  that  we  can  proceed  recursively 
with  ft  — ►  oo  unless  we  have  ft  —  oo  for  some  v  —  P°  <  oo.  The  argument  also 
shows  that  [ xn ,  jn]  — ►  [x°,  j°]  on  [0,  T°]  for  each  v  (up  to  v  —  v  if  v  <  oo  ) 
where,  of  course,  this  refers  to  the  resulting  subsequence  remaining  after  all  the 
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extractions  to  this  point.  Clearly  this  convergence  carries  over  to  the  subsequence  (of 
the  orginal  sequence)  obtained  by  a  Cantorial  diagonal  construction  since,  from  some 
//-dependent  point  on,  this  coincides  with  the  subsequence  directly  associated  with 
arriving  at  t°u.  The  definition  of  convergence  on  [0,  oo)  shows  that  we  must  necessarily 
have  jn(m)]  — ►  [xQ,j°\  on  [0,  oo)  for  the  diagonal  subsequence,  as  asserted  and, 

by  Theorem  1,  we  then  also  know  that  [x°,  j°]  is  a  solution  of  E°.  □ 

We  turn  next  to  an  existence  theorem  for  global  solutions,  i.e.,  on  [0,oo).  We 
first  provide  a  cautionary  example,  showing  that  even  in  a  case  with  J  —  {1,2}  with 
each  7 Tj  asymptotically  stable  it  is  (surprisingly?)  possible  for  the  switching  system 
to  support  ‘blowup’  —  solutions  which  ‘escape  to  infinity’  in  finite  time. 

Example  3:  Take  X  —  5?2  and  let  7Ti,7T2  be  the  solution  operators  associated  with 
the  ordinary  differential  equations: 


(4.7) 


(i)  i  =  (l  +  l*|2) 


10 

-4 


(ii)  x  =  (1  +  |z|2) 


1 

5 


5 


respectively.  Clearly  the  solution  paths  in  3?2  are  the  same  as  for  the  linear  equations: 

(i)  ii  —  u  +  5v  +  10, 

,  .  v  =  —  5u  —  2v  —  4; 

(ii)  u  =  u  —  5v  +  10, 

v  =  5u  —  2v  +  4 


which  are  easily  seen  to  be  exponentially  stable:  the  characteristic  exponents  for  (  4.8) 
are  [—1  ±  4*]/2.  For  any  starting  point,  then,  (  4.7)  has  bounded  solutions  so  the 
velocities  are  also  bounded  and  the  solutions  go  exponentially  to  the  same  attractors: 
(0,-2), (0,2). 

Now  take 

Ri  :=  {(u,v)  :  v  <  — 1},£2  :=  {(u,t>)  :  v  >  1} 

Si, 2  :=  {(iz,-l)},S2>i  :=  {(u,l)}- 

One  easily  sees  that  every  solution  of  the  resulting  switching  system  E  will  alternate 
modes  with  switching  points  alternately  on  v  —  ±1  so  (eventually)  one  keeps  |i>|  <  1. 
If  one  then  has  u(t)  >  0,  we  subsequently  have  u  >  5(1  +  u2),  giving  blowup  before 
t  +  7t/20.  (Similarly,  u(t)  <  —16  gives  w  <  —  (1  +  u2).  ) 

Note  that  if  we  had  used  the  linear  equations  (  4.8),  instead  of  (  4.7),  to  define  the 
modes  7Ti,  7t2  then  we  could  not  get  blowup  in  finite  time  but  nevertheless  would  obtain 
instability  (unbounded  solutions)  for  a  switching  system  comprised  of  exponentially 
stable  modes.  □ 


16 


Theorem  3:  Let  E  :=  [J ,  {717},  {Sy,*}]  be  a  switching  system.  Suppose  the  index 
set  J  is  finite  and  each  Sj,k{j  7^  k)  is  compact.  Then  for  any  consistent  initial 
data  (£,j)  £  X  x  J  (i.e.,  £  Rj  )  there  is  global  regular  solution15  [x(-),j(-)]  with 

®(0)  =  im  =  J- 

PROOF  :  Simply  follow  the  condition  (  3.4)  (ii)  in  developing  the  evolution  of  x(-) 
with  switching  as  permitted/required  by  (  3.6),  choosing  almost  arbitrarily  when 
nonunique  choices  may  occur.  We  show  that  switching  times  will  ‘almost  automat¬ 
ically’  be  isolated  and  that  the  construction  cannot  stop  (i.e.,  one  can  neither  have 
finite  escape  time  nor  an  impasse  for  which  no  admissible  continuation  exists). 

We  set  £0  :=  £,Jo  =  j(0)  j.  This  is  permissible  since,  by  assumption,  £  0  Ry. 
It  is  possible  that  £  £  for  some  k  ^  j,  in  which  case  we  can  choose  to  switch 
immediately,  making  0  =  to  a  switching  time  and  ji  =  k  ^  jo  =  J.  If  one  had 
£  £  dZj  with  TTj  entering  Rj  at  £  then  an  immediate  switch  is  mandatory  by  (  3.6) 
(ii)  and  some  switch  is  permissible  by  (  3.3)  (iii);  else  one  could  permissibly  choose 
ji  =  Jo  =  3- 

Now  proceed  in  the  mode  7Tj(j  =  j\)  until  at  time  t  one  arrives  at  some  Sjtk{k  j). 
Either  £*  :=  x(F)  7Tj{t  —  to,  £0)  is  in  dRj  and  ir j  enters  Rj  at  £x  so,  as  above,  a 
switch  is  mandatory  by  (  3.6)  (ii)  and  available  by  (  3.3)  (iii)  or  an  optional  choice 
is  available  and  one  can  choose  either  to  continue  in  tv j  or  to  switch,  making  t  a 
switching  time  t%. 

Next,  we  observe  that  there  is  an  absolute  minimum  interswitching  time  f  >  0 
before  mandatory  switching16.  To  see  this,  note  that  if  we  make  a  switch:  j  ^  k  it 
must  be  at  some  £  £  Sj:k •  Then  the  evolution  proceeds  in  mode  7Tj t  and  a  subsequent 
switch  cannot  be  mandatory  unless/until  nk(t,^)  =:  £  would  be  in  dRk  for  some  t 
(and  7T k  enters  Rk  at  £  ).  By  (  3.3)  (i),  the  closed  sets  Sjtk  and  Rk  are  disjoint  so,  by 
the  continuity  of  TTk,  there  exists  r  =  v* ( <£)  >  0  and  M  =  .V(£)  such  that  Trk(t,  £)  £  Rk 
for  (t,£)  £  [0,r]  x  j/.  By  the  assumed  compactness  of  Sjtk,  we  can  cover  Sjtk  by  a 
finite  number  of  these  neighborhoods  and  let  Tjtk  be  the  smallest  of  the  associated  r's. 
The  minimum  of  these  (over  the  finite  set  {(j,  k)  £  J  x  J  :  j  ^  k}  )  is  then  f. 

The  assumptions  made  for  this  theorem  do  not  imply  disjointness  of  Sj,k  and 
Sk,i  ( k  7^  i,j),  etc.,  so  it  would  be  conceivable  that  the  geometry  would  permit  such 
voluntary  switching  as  to  violate  the  condition  (  3.4)  (i).  This  is  the  point  of  the  use  of 
“almost”  in  the  first  paragraph  of  this  proof:  we  must17  restrain  our  voluntary  choices 
so  as  to  avoid  violation  of  (  3.4)  (i).  Certainly  the  minimum  time  f  for  mandatory 

15Indeed,  we  see  that  (almost)  every  solution  constructed  following  the  switching  rules  (  3.6)  is  regular 
(satisfying  (  3.4)  (i))  and  global. 

16Note  that  in  many  contexts  it  is  reasonable  to  have:  Sj,k  C  dRj  for  k  ^  j.  In  this  case  one  would 
have  the  same  minimum  f  for  voluntary  switching  as  is  obtained  here  for  mandatory  switching. 

■^Alternatively,  we  could  work  with  a  more  general  notion  of  solution  but  this  leads  to  technical 
complications  and  we  do  not  discuss  this  here;  see  Section  6,  however. 
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switching  means  that  it  is  always  possible  to  avoid  violation  of  (  3.4)  (i)  —  e.g.,  one 
could  choose  to  switch  only  when  this  is  mandatory. 

Clearly  the  construction  of  [x(-), ./(•)]  then  proceeds  for  t  -*■  oo.  Note  that  the 
number  of  switching  times  may  be  infinite,  with  x(-)  defined  inductively  on  the  inter¬ 
switching  intervals  by  the  appropriate  mode:  x(t)  =  7Ty(i— t„-i,  £„_i)  for  t„_  x  <t  <tv 
with  £„  then  given  as  x(tv).  On  the  other  hand,  there  may  be  a  last  switching  time  t 
with  x[t)  —  —  £)  for  t  <  t  <  oo  (  assuming  this  never  would  give  x(t)  £  trig¬ 

gering  a  mandatory  switch).  In  either  case,  this  defines  [rn(-) , j (•)]  :  3i+  ->  X  X  J  and 
it  is  clear  from  the  construction  that  (assuming  such  restraint  in  voluntary  switching 
as  may  be  needed  to  avoid  violating  (  3.4)  (i))  any  such  [x,j]  is,  indeed,  a  solution  of 
E  as  desired.  We  may  note,  also,  that  all  solutions  are  obtainable  in  this  way.  □ 

One  might  feel  that  the  restriction  to  finite  J  and  compact  Sj,k  for  j  ^  k  avoids 
the  difficulties  observed  in  connection  with  Example  3.  Indeed,  Theorem  2  shows 
that  ‘escape  to  infinity  in  finite  time’  is  then  impossible.  It  is  tempting  to  conjecture 
that,  with  these  restrictions,  if  each  component  mode  has  the  boundedness  property : 

(4.9)  for  each  solution  the  set  (x(f)  :  t  £  9?+}  is  bounded 

then  the  same  would  be  true  for  the  switching  system.  We  will  see  that  this  need 
not  be  true  even  if  the  state  space  X  is  taken  to  be  locally  compact:  we  provide  a 
counterexample  with  X  —  dt2  and  two  modes. 

Example  4:  Begin  by  defining  a  mode  7Ti  as  the  solution  operator  for  a  differential 
equation  :  x  =  f(x)  for  x  =  ( u,v )  £  3 J2  =:  X.  The  direction  of  (the  2-vector)  f(x ) 
will  be  specified  by  describing  the  integral  curves;  the  speed  |/(x)|  along  the  curves 
can  be  specified  independently.  We  first  set 

S  :=  {u  <  -2}  U  {|u|  >  1}  U  {(u,0)  :  |u|  <  2}. 

This  will  be  the  set  of  stationary  points  for  the  flow  (/(z)  =  0)  and  elsewhere  we  set 

|/(x)|  :=  dist  (x,  S). 

The  direction  of  /  is  irrelevant  (undefined)  on  S,  particularly  for  { u  <  —2}  and 
{|u|  >  1},  and  need  only  be  determined  on 

{(u,v)  :  u  >  —2,  |v|  <  1}  \  {(w,0)  :  u  >  —2}. 

For  u  <  2  we  take  f(u,v )  to  be  of  the  form:  (a,0)  with  sgna  —  —sgnv.  Finally, 
we  fill  out  the  half-strip  {(u,u)  :  u  >  2;  |v|  <  1}  with  the  family  of  integral 
(C(A)  :  0  <  A  <  1}  given  by 

C(A)  :=  {(2  +  [V\  —  v][A  +  v\/\l  —  A],  v)  :  —A  <  v  <  V\}. 


curves 
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Each  point  x  =  ( u,v )  in  the  (open)  half-strip  is  on  a  unique  integral  curve  C(A)  and 
we  take  the  direction  of  f(x)  to  be  tangent  to  that  curve,  oriented  so  that  motion 
along  the  curve  is  counterclockwise;  see  Figure  3.  This  defines  7Ti. 


Now  define  another  mode  7r2  by  setting  7T2(t,  Q  :=  — 7Ti(t,  — £)  on  5R+  x  SR2.  Clearly 
we  have  (  3.1)  for  each  7r  =  7ry  (j  E  J  :=  {l,2})and  it  is  easy  to  see  that  (  4.9)  holds 
for  each  mode  since,  for  7Ti,  the  set  {x(t)  :  t  E  S?+}  for  a  solution  consists  either  of  a 
single  point  in  S  or  of  (part  of)  a  segment  {(u,\/A)  :  — 2  <  u  <  0}  with  0  <  A  <  1 
together  with  (possibly)  (part  of)  the  curve  C(A)  together  with  (possibly)  (part  of) 
the  segment  {(u,  —A)  :  —2  <  u  <  0}.  We  set 

£i  ~  {(u,v)  :  -2  <  u  <  — 1;  0  <  v  <  1}  ,  Si)2  :=  dUi  ; 

R2  •—  {(«,v)  :  1  <  u  <  2,  -1  <  v  <  0}  ,  52,i  :=  dZ2 
to  complete  the  specification  of  E. 

Now  consider  the  (unique)  solution  of  £  starting  in  mode  7Ti  (i.e.,  j(0)  =  j0  =  j\  = 
1  )  at  a:(0)  =  £0  =  (0,  —  Aj)  with  0  <  Ai  <  1.  The  state  moves  (in  mode  )  to  the 
right  until  reaching  (2,— Aj),  then  lo.ops  around  on  C(Ai)  until  reaching  (2,A2)  with 
A2  :=  \/Ai,  then  moves  left  until  reaching  the  switching  point  :=  (— 1,A2).  All  of 
this  takes  finite  time  so  the  speed  is  bounded  away  from  0  (since  the  path  is  bounded 
away  from  S).  Now  the  state  continues  moving  left  (in  mode  7r2)  until  reaching 
(— 2,A2),  then  follows  —  C(A2)  around  until  reaching  (— 2,— A3)  with  A3  :=  then 
moves  right  until  reaching  the  switching  point  f2  :=  (1,— A3).  Etc.  The  complete 
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path  of  this  solution  is  an  expanding  ‘spiral’  which  is  composed,  alternately,  of  curves 
C(A2n-i)  and  —  C(A2n),  connected  by  straight  segments  along  v  =  ±A„.  Note  that 
A„+i  =  [v  =  1,2,...)  so  Aj,  — >•  1.  Since  each  curve  C(Ai/)  includes  the  point 
(2  +  A^/2/[l  —  A„],  0),  this  shows  that  the  complete  path  is  unbounded:  (  4.9)  fails  for 
E.  □ 

5.  Linear  Switching  Systems 

As  a  particularly  interesting  and  important  case,  we  consider  certain  bimodal 
systems  (J  =  {1,2})  with  linear  dynamics. 

Let  X  be,  e.g.,  a  real  Banach  space  and  suppose  A  is  the  infinitesimal  generator 
of  a  Co  semigroup  S  of  linear  operators  on  X: 

(i)  S(t  +  a)  =  S(t)S(s),  S(0)  =  1; 

(5.1)  (ii)  t  i — *  S (t)£  continuous  in  t  for  each  £  6  X  ; 

(iii)  ||  S  (4)  ||  <Mewt  forte  $R+. 

We  now  assume  the  two  modes  are  given  by 

(5.2)  x  =  Ax  +  uj  ( j  =  1, 2) 

where  Ui,u2  are  specified  (constant)  elements  of  X  so  the  standard  variation  of  pa¬ 
rameters  formula  gives 

(5.3)  iv j (t ,  £)  =  S(t)£  +  S (s)ds]uj 

We  also  suppose  A  €  X*  is  specified  (with  A  ^  0)  and  set 

(5.4)  Sj  :={£€*:  (A,  <  «x}, 

£2  :=  {£  £  X  :  (A,  £)  >  a2}, 

Si, 2  •—  {(A,  £}  =  oci}  =  dZi  ,  52,i  =  {(A,  £)  =  a2}  =  dR2 

where  <  a2  are  specified  scalar  values.  A  switching  system  E  specified  as  in  (  5.3, 
5.4)  will  be  called  a  linear  switching  system.  What  characterizes  linear  switching 
systems  is  the  representation  for  solutions: 

(5.5)  x(t)  —  S (t)£  +  f  S (t  =  s)uj(s)ds 

J  o 

=  Xo(t)  +  [  cr1(s)[S(t  —  s)«i]ds  +  [  a2(s)[S(t  —  s)u2]ds 
Jo  Jo 

where  x0(i)  :=  S (t)£  with  £  :=  x(0)  and 

°k{t)  :=  {1  if  j{t)  —  k ;  else  0}  ( k  =  1,2) 
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Although  the  switching  surfaces  SJtk  are  not  compact,  as  assumed  for  Theorem  3,  we 
see  that  one  still  has  global  existence: 

Theorem  4:  Let  E  be  a  linear  switching  system.  Then  for  any  consistent  initial 
data  (£,  j)  the  set  of  global  solutions  is  non-empty  and  equicontinuous  with  a  lower 
bound  r(T')  on  lengths  of  interswitching  times. 

PROOF  :  The  argument  for  existence  is  exactly  as  in  the  proof  of  Theorem  3  except 
for  the  discussion  of  the  minimum  for  interswitching  intervals.  Since  the  switching 
surfaces  Sij2,  £2,1  are  parallel  hyperplanes,  we  see  that  any  interswitching  interval 
must  permit  time  for  a  transit  across  the  gap. 

From  the  representation  (  5.5)  we  can  see  that  any  such  function  (  i.e.,  whether 
or  not  one  obeys  the  switching  rules  (  3.6)  )  satisfies 

(5.6)  \x(t)  -  x(t)|  <  \x0[t)  -  x0{t)\  +  f  \S(t  -  s)uj^)\ds 

J  t 

<  |x0(t)  -  x0(t)|  +M\t  -  t\ 

where 

M  =  M(T')  :=  Mmax  {[ttjj,  |u2|}raax  {l,ewT’}. 

If  (t,t)  is  an  interswitching  interval,  then  (  5.4)  gives 

a2  -  ai  =  j  {A ,x(t)  -  x(t))  ]  <  |A|  jx(t)  -  x(t)\. 

Given  any  V  <  00,  x0(-)  is  uniformly  continuous  on  [0,  T']  so  there  is  some  r0  =  ro(T') 
such  that 

| x0(t)  -  x0(t)|  <  (/?  -  a)/ 2  jA|  if  (t-t)  <  t0 
for  0  <  t  <  t  <  T';  we  aslo  set 

n  =  ri(T')  :=  (a2  ~  ai)/w  |A|  M. 

Then  (  5.6)  gives  the  lower  bound 

(5.7)  (f  —  t)  >  min{r0,Ti }  =:  r(T') 
as  well  as  the  equicontinuity.  □ 

Remark  2:  Suppose  ^?(s)  is  a  nondecreasing  bound  for  |S(s)£  —  £|  so  (  5.1)  (ii)  gives 
<£>(.s)  — *  0  as  s  — >  0;  let  ip  be  the  inverse  function  of  <p.  Since 

\xQ{t)  -  x0(OI  <  ||S(f)||  |S(f-  t)  i-i\<  Mewt(p{t  —  t), 

we  can  use  ^([0:2  —  aq]/2  |A|  Mmax{l,  ewT'})  as  tq{T').  In  particular,  if  £  6  P(A)  we 
can  take  £>(s)  :=  MjA^|s  for  small  s  (any  M  bounding  Me“8  )  and  this  gives  r(-) 
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decaying  at  worst  like  e  U>T>  if  u>  >  0  and  fixed  (independent  of  T1)  in  the  stable  case: 
uj  <  0.  Since  what  we  are  really  estimating  is 

<  A,  xo(t)  —  £o(t)  >=<  S *(t  —  t)A  —  A,  S (t)£  >, 

we  could  proceed  rather  similarly  if  we  bounded  \S*(t  —  t)X  —  A|  instead.  In  particular, 
if  A*  generates  an  adjoint  semigroup  S*(-)  on  X*  and  if  A  £  D(A*),  then  we  obtain 
the  same  ‘at  worst  exponential’  decay  rate  for  r(-),  now  for  arbitrary  £  £  X.  □ 

We  call  a  linear  switching  system  £  stable  if  the  defining  semigroup  S(-)  is  (expo¬ 
nentially)  stable:  oj  =  —a  <  0  in  (5.1)  (iii).  Note  that  in  this  case 

|£|*  :=  sup{eat\S{t)ti\  :  t  £  9i+} 

defines  an  equivalent  norm18  on  Af(|£|  <  |£|*  <  M|£|)  and  use  of  this  norm  makes 
M  =  1  in  (  5.1)  (iii).  Without  loss  of  generality  we  henceforth  assume  that  M  =  1  so 

||S(t)||  <  e~at 

with  a  =  —u  >  0.  If,  in  addition,  one  has 

(5.8)  S(er)  is  compact  [  so  S(t)  is  compact  for  t  >  a] 

for  some  o  >  0,  then  we  call  £  compact  but  we  note  that  stability  alone,  without 
(  5.8),  already  gives  certain  compactness. 

Lemma  1:  Let  S(-)  be  stable  and,  for  U  C  X,  set 

So  —  So (U)  :=  {0}  U  {S(t)u  :  t  €  9i+,u  G  U}, 

It  =  It{U)  :=  {J  S (t)u(t)dt  :  u(-)  measurable,  U  —  valued}. 

Suppose  U  is  compact.  Then  So  is  compact  and,  for  each  T  >  0  (including  T  =  oo  ), 
It(U)  is  precompact  in  X. 

PROOF  :  Note  that  (  5.1)  gives  continuity  of  S [t)u  jointly  in  [t,u)  so  the  image  of 
[0 ,T]  x  U  is  compact  for  any  T  <  oo.  On  the  other  hand,  the  tail  is  covered  by  any 
e  -ball  at  0  if  T  >  [log  p/e\/a  where  p  :=  max{\u\  :  u  £  11}  <  oo.  This  makes  So 
precompact  and  one  easily  sees  that  it  is  closed  in  X,  hence  compact. 

^Alternatively,  we  could  take 

/»  oo 

let  .:=\J  e2a‘|S(t)e|  2dt}V* 

to  retain  a  Hilbert  space  structure  for  X  if  it  is  originally  Hilbertian. 
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For  It  we  consider  T  <  oo  first.  Then  one  easily  sees  that  It(U)  C  T$(U),  hence 
is  precompact.  For  T  =  oo  we  show  loo  totally  bounded.  For  any  e  >  0  we  cover  IT> 
by  e/2-  balls,  taking  T'  >  [log  (2 fi/ea)]/a  so,  uniformly, 

r  oo  /*oo 

|  /  S(t)u(t)dt|  <  /  e~at [idt  <  e/2 

Jr'  Jr' 

Then  Joo  is  contained  in  the  e-ball  cover  with  the  same  centers:  each  integral  in  I ^ 
is  the  sum  of  one  in  Iy»  and  a  small  tail.  □ 

Theorem  5:  Let  E  be  a  stable  linear  switching  system.  Then  there  is  a  bounded 
invariant  set  B  C  X  such  that  ( the  state  component  of )  every  solution  of  E  eventually 
enters  and  stays  in  B.  For  any  compact  set  U  C  X  there  is  a  compact  convex  set  S  (U) 
such  that  (the  state  component  of)  every  solution  which  starts  in  U  will  lie  entirely 
in  S(U).  Finally,  if  E  is  compact  then  B,  above,  can  be  taken  to  be  compact. 

PROOF  :  From  (  5.5),  any  solution  of  E  has 


r  oo 

x(t )  =  S(t)x(0)  +  /  S  (s)u(s)ds 
Jo 


with 


u(s)  :=  for  s  <  t;  0  for  s  >  t}. 


Thus,  x(t)  will  be  in  S(U)  for  all  such  x(-)  with  x(0)  6  U  and  for  each  t  &  9i+  if  we 
set 


S(U)  :=  co[S0(W)  +  I*] 


where  I*  :=  /Oo({tti5«2?0})-  Both  So{U)  and  I*  are  compact  by  Lemma  (Ll)  so  this 
S(U)  is  compact  and  convex  as  desired. 


The  same  analysis  gives  x(t)  €  [Be  +  I*]{Be  :=  {£  G  X  :  |£|  <  e})  when  t  > 
[log  (\x(0)\l e)\l a  so,  e.g.,  [Si  +  I*\  B  is  globally  attractive.  This  B  will  not,  in 
general,  be  invariant  so  we  take 


(5.9)  S  :=  {S(t)|  +  Jt({tti,u2})  :  I  G  B,t  G  8J+}. 

Taking  t  =  0,  we  see  that  S  C  S  so  S  is  globally  attractive;  S  is  obviously  bounded 
and  we  need  only  show  invariance.  For  £  G  S  we  have,  from  the  definition  (  5.9),  that 


£  =  S(*)f  +  f  S(i  -  s)u3^ds 

J  0 


for  some  £  G  B,t  £  3?+,  measurable  ,;(•)  :  [0,t]  — >  J.  Any  solution  [x(-),j(-)]  with 
x(0)  =  £  will  have,  by  (  5.5), 

x(r)  =  S(r)[S(t)£  +  f  S(t  -  s)uS(3)ds}+  f  S(r  -  e)uj{e)dg 

J  0  •'0 

a  rt+r 

=  S(t  +  r)£  +  /  S(t  +  r  -  s)v,j(s)ds 

J  o 


where  we  have  set 


i(s)  := 


j(s)  on[0,<], 
j[s-t)  on  [t,t  +  r]. 


This  shows  x(r)  6  B  for  any  such  z(-),  any  r  E  3f?+ ,  i.e.,  B  is  invariant. 

A 

Finally,  suppose  E  is  compact,  i.e.,  (  5.8).  Clearly,  if  x(t )  £  S  for  any  solution 
[*(•), i(0]  an<^  some  t  €  3i+,  then  x(t  +  cr )  is  in  T0  :=  [S(er)B  +  1*]  which  is  precompact 
by  (  5.8)  and  Lemma  1;  hence,  the  compact  convex  set  T2  :=  coTi  is  a  global  attractor. 
Again,  this  may  not  be  invariant  but  we  can  introduce 


Ts  :=  {S(*)£  +  It{co{uu  u2})  :  |  G  T a,  f  €  3 ?+} 

which,  as  for  B  above,  is  an  invariant  global  attractor.  It  is  not  difficult  to  verify 
from  its  form,  noting  the  definition  of  r2,  that  r3  is  convex.  Since  T3  C  [So(T2)  +  I*], 
it  is  precompact.  We  set  T  :=  T3  and  note  that  T  is  obviously  convex,  compact,  and 
globally  attractive.  To  see  its  invariance,  note  that  if  [:r(*),  ./(•)]  is  any  solution  of  E 
starting  at  £  £  T  (so  one  has  £k  — >  £  with  £k  £  T3  ),  then  for  any  f  G  9J+ 

xk(t)  :=  S(t)£k  +  f  S(t  -  s)uj(a)ds 

J  0 

->  S(t)t+[  S (t  -  s)uj^ds  =  x(t) 

J  o 

by  the  continuity  of  S  (t).  On  the  other  hand,  the  proof  of  invariance  of  T3  shows 
%k{t)  £  r3  so  x(t)  6  f 3  =  r.  □ 

We  conclude  this  section  by  considering  solutions  on  all  of  9J  for  a  stable  linear 
switching  system  E.  For  such  a  solution  [z(-), ./(■)]  :  3?  — >  X  x  J,  we  refer  to  the 
restriction  of  j’(-)  to  (— oo,t\  as  the  switching  history  (at  time  t)  of  the  solution;  it  is 
convenient  to  represent  this  by  a  {0,  l}-valued  function  77(f)  on  3?+: 


(5.10)  fr(i)](s)  =  v{t’,s)  :=  j(t  -  s)  -  1. 

At  least  for  the  class  of  bounded  solutions  (i.e.,  with  x(-)  £  L°°(5i  — »■  X)),  we  will 
see  that  the  restriction  of  x(-)  to  (— 00,  t]  can  be  recovered  from  This  permits 
a  canonical  representation  in  terms  of  a  system  as  in  Example  2;  abstractly,  the 
dynamics  and  switching  are  characterized  by  the  function  appearing  in  (  3.13)  and 
the  values  of  Ci,C2  in  (  3.17). 

Theorem  6:  Let  E  be  a  stable  linear  switching  system  as  in  (  5.3),  (  5.4)  such  that 
each  Rj  is  globally  attractive  for  7 ly.  Then,  if  [»(•), j(-)]  is  a  bounded  solution  of  E 
on  (— 00,  f]  and  77 (f;  s)  is  given  by  (  5.10),  we  have 
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r  oo 

with  Vj  :=  /  S (s)ujds. 

Jo 

Next,  let  E°  be,  as  in  Example  2,  defined  by  (  3.10),  (  3.17)  with 
(5-12)  ip(t)  :=  c  <  A, S(t)[u2  —  Mi]  >  t  € 

with  c  :=  1/  <  A,v2  —  Ui  >  and  using  suitable  constants  Ci,C2  in  (  3.17).  Let 
[rc(-),y(-)]  :  8?  — ►  X  x  J  be  a  bounded  solution  on  (— oo,0]  =:  9i~  of  E;  define 
?/(•)  :  3?+  — >  X°  by  y{t)  r](t;-)  as  in  (  5.10)  and  let  j  be  the  restriction  of  j(-)  to 
9J+.  Then  [x(-), .?(•)]  is  a  solution  on  of  E°  and  a:(-)  is  given  on  3i+  by  (  5.11). 

PROOF  :  Define  x(t)  on  (— oo,f|  by  the  right  hand  side  of  (  5.11)  with  t  (variable) 
replacing  t.  Using  (  5.5)  we  have 

x(t)  =  S(t  —  t)x(t)  +  ^  S (t  —  s)uj(3)ds 
for  t  <  t  as  the  restriction  of  [m(-),  j  (•)]  is  a  solution  on  [t,t\  of  E.  Since 
Uj(s)  =  «1  +  [j{s)  -  l][ua  -  tti]  =  ux  +  —  s)[u2  -  «i], 

a  bit  of  manipulation  gives 

x(t)  —  x(t)  =  S  (t  —  t)[x(t)  —  x(t)} 

for  arbitrary  t  <  t.  As  x(-)  is  bounded  by  assumption  and  one  easily  verifies  that  x(-) 
is  also  bounded,  this  gives 


|z(t)  —  x(t)|  <  e 

for  some  M  and  arbitrarily  large  [t  —  t),  whence  x(t)  =  x{t)  as  desired. 

Our  choice  of  c  for  (  5.12)  gives  the  normalization  (  3.12)  (i);  if  we  take  (p(t)  := 
c|A|e~“*  in  (3.8),  then  (  3.9)  and  (  3.12)  (ii)  hold.  The  condition  (  3.12)  (iii)  can  be 
obtained  along  the  lines  of  Remark  2  if  either  [w2  —  ui]  <E  V(A)  or  A  €  D( A*)  but  we 
always  have 

(5.13)  1^(0  <^(«)7(*-'S) 

where  7  :  SR+  — >•  3?+  with  ^(r)  — >  0  as  r  — >  0, 

taking  <p  as  above  and  ■j(r)  :=  j  [S  (r)  —  l][u2  —  ui]|.  Going  over  Example  2  we  see 
that  (  3.12)  (iii)  was  used  only  to  obtain  a  uniform  lower  bound  on  the  interswitching 

A 

intervals  through  Lipschitz  continuity  of  0.  We  would  now  have  (  3.16)  replaced  by  a 
uniform  continuity  estimate 
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and  this  again  gives  tv+ 1  —  tv  >  Tmin  as  in  (  3.18).  We  will  be  taking 

(5.14)  Cj  :=  c[«j—  <  A,^  >]  (j  =  1,2) 

for  (  3.17).  One  easily  sees  that  Vj  is  globally  attractive  for  7 r,-  (j  =  1,2)  so  the 
assumption  that  Rj,  as  given  by  (  5.4),  is  attractive  for  7 Xj  implies,  with  our  definitions, 
that  0  <  Ci  <  C2  <  0.  We  see  from  (  5.11)  that 

roo 

<  A ,x(t)  >  =  <  A,vi  >+  r](t;s)ip(s)ds/c, 

Jo 

a  r  o° 

0(t)  —  /  r)(t]  s)ip(s)ds  =  c[<  A ,x(t)  >  —  <  A,iq  >] 

Jo 

so  (  5.14)  gives 

(5.15)  6{t)  —  Cj  <=$■<  A,x(t)  >=  a.j 

with  y(t)  =  entering  the  region  under  7r°  (i.e.,  in  the  senses  of  (  3.17), 

(  3.10))  precisely  if  x[t )  were  entering  Rj  under  tt j  (in  the  senses  of  (  5.3),  (  5.4)). 
Thus  the  switching  rules  for  E  and  E°  exactly  correspond,  specifying  precisely  the  same 
permissible  and  mandatory  transitions.  The  validity  of  (  5.11)  and  the  dynamics  for 
the  switching  components  then  give  the  desired  correspondence19  of  solutions.  □ 

6.  Comments  and  Discussion 

It  was  already  mentioned  above  that  more  detailed  consideration  of  certain  direc¬ 
tions  for  further  investigation  will  be  presented  in  [6],  [7].  We  note  here,  quite  briefly, 
some  other  possible  extensions  of  the  notions  presented  here  and  possible  areas  for 
further  investigation. 

The  most  obvious  of  these  would  be  consideration  of  time- dependent  problems. 
This  could  involve  either  time-dependent  modes  or  time-dependent  sets  {Sy,*}  or 
both.  As  noted,  the  usual  trick  of  absorbing  the  time-dependence  into  an  augmented 
state  handles  the  local  theory  but  unduly  complicates  the  consideration  of  global 
(asymptotic)  questions.  No  essentially  new  ideas  seem  to  be  involved  here. 

Slightly  more  delicate  would  be  the  treatment  of  solutions  for  which  one  admits 
the  possibility  of  a  limit  of  switching  times,  a  possibility  we  have  eschewed  here  in 
restricting  consideration  to  ‘regular  solutions’.  If  t  would  be  a  limit:  tv  — >  t—  with 
each  t„  an  isolated  switching  time  for  a  switch:  jv  ^  consider  the  set  of  pairs 

P  :=  nP{(i,i+i)  :  v  >  v} 

19Note  that  E  may  have  unbounded  solutions  on  9J  for  which  (  5.11)  does  not  apply  and,  of  course, 
E°  has  solutions  on  SR+  which  do  not  correspond  to  any  solution  on  SR+  of  E.  The  correspondence 
is:  {  solutions  of  E°  on  3?}  <-►  {  bounded  solutions  of  E  on  3?  }. 
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and  note  that  this  must  be  nonempty  if  we  have  (  3.3)  (iv).  Then,  without  attempt¬ 
ing  to  assign  a  value  to  j(t),  it  would  be  plausible  to  accept  a  solution  with  j(-)  =  k 
on  an  interval  (f,  ?)  for  any  k  such  that  some  (j,  k)  is  in  P  since  j(-)  must  take  the 
value  j  infinitely  often  as  tv  — ►  t—  and  x(t)  must  be  in  Similar  considerations 

would  apply  to  consideration  of  j(-)  constant  on  an  interval  (?,i)  and  then  switching 
infinitely  often  to  the  right  of  t.  These,  however,  would  be  only  the  simplest  possibil¬ 
ities  and  it  is  not  entirely  clear  how  to  define  ‘solution’  so  as  to  omit  or  significantly 
weaken  the  condition  (  3.4)  (i). 

An  interesting  possibility,  also,  would  be  to  weaken  the  implicit  assumption  that 
X  is  specified20  before  introducing  the  modes  {7 r,-  :  j  £  J}.  One  could,  alternatively, 
associate  a  state  space  Xj  with  each  tt j  and  then  ‘glue  these  together  at  the  switching 
sets  Sjtk,  i.e.,  one  would  have  (continuous)  functions  Tjtk  •  Xj  Z>  Sjtk  — ►  Xk  and 
would  have  distinct  left-  and  right-hand  limits  for  x(-)  across  a  switch:  j  ^  k  so 
x{t—)  =  £_  €  Sj,k  C  Xj  and  x(t+)  =  £+  :=  Tjik(£-)  £  Xk.  In  (  3.4)  (ii)  we  would  use 
£_  for  x(t)  or  £+  for  x(s)  in  the  obvious  way.  Note  that  each  7 \j  would  only  be  defined 
on  Xj  (which  might  be  thought  of  as  manifold  with  boundary  dRj  )  and  ‘7 Tj  enters 
Rj  at  £'  would  mean,  simply,  that  £  £  dRj  :=  (bdry  Xj  )  is  not  a  possible  initial 
value  for  7 r,;  this  involves,  also,  a  modification  of  the  notion  (  3.1)  of  ‘mode’.  One 
can  generalize  the  notion  of  switching  system  by  not  requiring  each  tt j  to  be  defined 
on  all  of  3i+  x  X.  Let  us  weaken  (  3.1)  to: 

(i)  for  each  £  £  X  there  exists  an  interval  I  (£)  C  3J+  with 
£  1 — *  I(£)  upper  semicontinuous:  if  t  £  I(£n)  for  any 
£„  ->  £  in  X,  then  [0, fj  C  /(£); 

(ii)  n(t,  £)  £  X  is  defined  for  £  £  X,  t  £  I  (£) 
with  7r(0,  £)  =  £; 

(iii)  t  £  1  (7 r(s,  £))  <=>■  (t  +  s)  C  1  (£)  and  then 
tt (t  +  «,£)  =  7r(t,7r(s,  £)); 

(iv)  if  £n  — *■  £  £  X  and  t  £  1  (£„),  then 
tt(*,  £»)  ->  ?r(-,  £)  uniformly  on  [0 ,t\ 

Now  suppose  we  are  given  an  index  set  J,  a  family  of  state  spaces  {Xj  :  j  £  J} 
and  a  family  of  modes  {7 \j  :  j  £  J  each  acting  on  the  corresponding  Xj  as  in  (6.1) 
with  the  time  intervals  now  denoted  by  I/(£)  for  £  £  Xj,  j  £  J.  Next,  suppose  we 
are  given  sets  Sj,k  C  Xj  (possibly  empty)  for  k  ^  j  and  maps  Tjtk  :  Sjik  — *  Xk-  We 


(6.1) 


20 In  connection  with  this,  one  could  also  consider  j-dependent  index  sets:  Jj  :=  {k  :  one  has 
(nonempty)  sets  Sj,k  C  Xy}.  It  seems  likely  that  one  could  always  construct  a  ‘universal’  J  and  a 
universal  X  =  U{X,  :  j  €  J}  in  such  a  way  as  to  reduce  this  more  general  notion  to  the  original. 
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assume: 

(i)  each  Sy.jt  is  closed  in  Xj]  each  Tjtk  is  continuous  from 
Sj'k  to  Xk\£k] 

(6.2)  (ii)  for  each  £  G  Xj  there  is  a  neighborhood  M  of  £  for  which 
{k  G  J  :  [ Sj^  fl  ^  0}  is  finite; 

(Hi)  £j  :=  {£  G  Xi  :  MO  =  {0}}  C  Sj  :=  U {Sj>k  :  k  #  j}. 

The  specification 

S  :=  {^.{(Xj,^)  :  j  G  %/}, {(5y,*, 2y,fc)  :j,keJ,k j}}, 
as  above,  is  then  a  switching  system. 

It  is  now  convenient  to  define  a  solution  of  S  (on  an  interval  [0,T]  with21  T  <  oo) 
as  given  in  terms  of  a  finite22  partition 


0  =  t0  <  h  <  . . .  <  tp  =  T. 


We  assume  we  give,  for  u  =  1, . . . ,  P, 

(6.3)  j  —  G  J  and  £„(■)  :  \ju— *  Xj 

satisfying,  for  u  =  1, . . . ,  u  : 

%i/[t)  —  ^j(j  '  ly—ly  £v)  On  [tj/— l,t„] 
with  j  =  jv  and  £}  :=  Xv(tv-i)  G  Xj 
(implicitly  this  requires  [0,  —  i„_i]  C  MO); 

xv{tu)  =:  O  G  Sjik  and  0+i  -  Tjik(0) 
with  j  =  jv,  k  —  jv+i  for  1  <  v  <  V  (  and  for  v  =  0 
unless  ji  -  j0  and  £f  =  £0). 

Note  that  we  have  [£(•),  j(-)]  well-defined  on  [0,t]  except  at  {tv  :  u  =  0, . . . ,  P—1},  with 
j(-)  J-valued  and  a:(-)  taking  values  in  the  appropriate  state  space  for  each  subinterval 
of  the  partition.  It  is  then  sightly  awkward  to  define  x{tu )  but  we  have  the  limits 
X{tv~)  =  0-1  and  X[tv+)  =  O' 

It  is  not  difficult  to  see  that  any  switching  system  in  the  earlier  sense  becomes 
a  switching  system  in  the  sense  just  defined  on  taking  Xj  :=  X  \  Rj,  each  Tjtk  to 
be  the  (suitably  restricted)  identity  map,  and  each  mode  7r j  ‘maximally  defined’  on 
3ft+  X  Xj  to  give  (  3.1)  from  (  3.3).  The  requirement  that  Tj,k  not  take  values  in  £k  is 

21  We  now  define  a  solution  on  5R+  by  requiring  that  the  restrictions  be  solutions  on  [0,  T ]  for  arbitrarily 
large  T. 

22 This  subsumes  the  regularity  condition  (  3.4)  (i)  from  the  original  definition. 


(6.4) 


(i) 


(ii) 
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just  a  weaker  form  of  the  earlier  requirement  in  (3.3)  (i)  that  Sjik  H  Rk  =  0  since  we 
recognize  Sk  as  corresponding,  for  the  original  definition,  to  the  set  of  points  in  dRk 
at  which  -Kk  would  enter  Rk-  The  solution  set  for  the  re-interpreted  switching  system 
(new  definitions)  corresponds  precisely  to  the  original  solution  set. 

One  advantage  of  this  reformulation  of  the  notion  of  switching  system  is  to  make 
it  easy  to  define  an  extension  of  a  switching  system.  We  will  call  E'  :=  { J1, . . .}  an 
extension  of  E  and  write  E'  D  E  (equivalently,  we  call  E  a  restriction  of  S’)  if 
each  is  a  switching  system  (in  the  sense  of  our  reformulated  definition)  with 

J  C  J',Xj  closed  in  X]{j  €  J  C  J1),  Sjtk  C  Sjtk, 

and  if  7 ry,  Tjtk  are  the  appropriate  restrictions  of  7r '• ,  Tj  fc.  It  is  easy  to  see  that  every 
solution  of  E  is  also  a  solution  of  E'.  Conversely,  if  we  take  Sj,k  to  be  k  PI  Xj 
whenever  j,  k  €  J  C  J' ,  then  any  solution  of  E'  which  stays  in  {Xj  :  j  £  J}  will 
also  be  a  solution  of  E.  In  particular,  whenever  we  would  have  an  invariant  set  (as 
in  Theorem  3)  we  could  cut  down  the  state  space (s)  in  the  obvious  way  to  obtain  a 
restriction  of  E  -  indeed,  this  is  essentially  what  it  means  to  have  an  invariant  set. 

Finally,  we  note  the  interesting  possibility  of  considering  two-point  boundary  value 
problems  and/or  (optimal)  control  problems  in  the  context  of  switching  systems  as  a 
generalization  of  the  usual  problems  for  differential  equations. 
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